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The effect of the site potential in an imperfect layer is studied in a d-wave layer super- 
conductor on the basis of the electron correlation. The site potential binds electrons to the 
imperfect layer, and then, the superconductivity of the imperfect layer is independent of that 
of the bulk. We found that the superconducting transition temperature of the imperfect layer 
becomes higher than that of the bulk owing to the effect of the site potential. In this situation, 
the Fermi surface of the imperfect layer has a strong nesting feature leading to increases in both 
the antiferromagnetic spin fluctuation and the density of states near the Fermi level, which are 
in favor of d-wave pairing. 
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Recently, a curious phenomenon has been observed ex- 
perimentally in the eutectic ruthenium oxide Sr2Ru04- 
Ru^""' and the Ce compound Celrlns.^ In experiments^'^ 
on these two materials, the superconducting transition 
temperature Tc^p observed using zero resistance is higher 
than the bulk transition temperature Tc^buik obtained 
by both the specific heat and the Meissner effect. In 
the p-wave superconductor Sr2Ru04-Ru, Tc^p is twice 
7c. bulk — 1-5 K.^ In the d-wave superconductor Celrlus, 
Tc,p = 1.2 K is also higher than Tc,buik = 0.4 K.^ 

The existence of Ru lamellae has been observed^ in 
Sr2Ru04-Ru. Although Tc^buik should be obtained by the 
superconducting transition of the entire system, T^^p be- 
comes higher than Tc^buik in the case where the current is 
networked between interfaces of lamellae due to the inter- 
face superconducting state with a higher transition tem- 
perature. The interface superconducting state is bound 
to the interface and is separable from the bulk states. 
While the lamellar structures have not been observed in 
Celrlus,^ the two compounds have common properties 
as follows: One is that the two materials are unconven- 
tional superconductors with layer structures. The others 
are that both Ru d- and Ce /-electrons form the strongly 
correlated electron systems with cylindrical Fermi sur- 
faces, and the mechanisms of the superconductivity in 
the bulk originate from the electron correlation with the 
2-dimensional property.^' Therefore, it is expected that 
the interface superconductivity also exists in Celrlus, 
similarly to the case of Sr2Ru04-Ru. It is natural to 
theoretically study the mechanism of interface supercon- 
ductivity with Tc^p higher than Tc^buik on the basis of the 
effect of the electron correlation. 

In our previous study, we investigated the increase 
in the superconducting transition temperature of d-wave 
pairing at the surface on the basis of the effect of the 
electron correlation in the layer system. We found that 
the transition temperature at the surface becomes higher 
than Tc^buik owing to the effect of the site potential at 
the surface atoms. The effect of the site potential in- 



creases both the antiferromagnetic spin fluctuation and 
the density of states at the surface, which are in favor of 
d-wave pairing. From these results, we suppose that the 
superconducting transition temperature becomes higher 
even at the interfaces of the imperfect layers. In an ex- 
periment on Celrlns,^^ the effect of an impurity was re- 
ported as the mechanism behind Tc.p being higher than 
?c,buik- Therefore, we propose the existence of imperfect 
layers formed by the impurities or defects in Celrlus, 
and then, we investigate the effect of the site potential 
in an imperfect layer on the superconducting tempera- 
ture of the d-wave pairing in this study. Hereafter, we 
call the superconducting state in the imperfect layer the 
imperfect-layer superconductivity. We apply the Green 
function theory^^ of the surface superconductivity^^ to 
the imperfect-layer superconductivity. 

As shown in our model in Fig. 1, we consider a layer 
system with an imperfect layer, which is embedded par- 
allel to the ab plane. This imperfect layer is assumed to 
comprise atoms with the difference in site potential W. In 
Celrlns, the electronic property mainly originates from 
the /-electrons^' of the Ce atoms occupying vertices 
of the tetragonal lattice. The cylindrical Fermi surface 
along the fc^-axis* is denoted by the bulk dispersion: 

E{kx, ky, kz) = 2i(cos k^ + cos ky) + 21^ cos kz — ^. (1) 

Here, the nearest-neighbor transfer integral in the ab 
plane is t. The nearest-neighbor transfer integral per- 
pendicular to layers is tz, which is much smaller than t 
since the 2-dimensional property is strong. We set t = 1 
as the energy unit and tz = 0.2. We fix the chemical 
potential /i = —0.54 at T = 0.0005, and the value of n 
gives an electron filling of 0.384 close to that of Ce-115 
compounds.^' The space symmetry is conserved in the 
ab plane, although the translational symmetry is broken 
in the z-direction owing to the presence of the imper- 
fect layer. We use the mixed Bloch-Wannier representa- 
tion and Green function theory. Our representation is 
based on the Bloch representation in directions parallel 
to the imperfect layer. The number of layers is denoted 
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Fig. 1. Multilayer system with an imperfect layer. The imperfect 
layer is positioned horizontally {ab plane) in the bulk layers. The 
gray and white circles indicate atoms of the imperfect layer and 
those of the bulk layers, respectively. The nearest-neighbor trans- 
fer integral in the ab plane is f. The nearest-neighbor transfer 
integral perpendicular to layers is t^. The difference in transfer 
between the imperfect layer and the bulk layers is defined by V. 
The value of the site potential in the imperfect layer is denoted 
by W. 



by m, and the imperfect layer is numbered by m = 0. 
Here, tz between the imperfect layer and the bulk layers 
at TO = ±1 is reduced by the difference V in transfer, 
which is fixed a,s V = — 0.25tz. The Green functions in 
this system including the imperfect layer are obtained 
using the Dyson-Gorkov equation: 

+G,f,_i(/c|| , u;i)V + G^+i (fc|| , , ^i) 

where cui = TrT{2l + 1) is the fcrmion Matsubara fre- 
quency. The momentum in directions parallel to the im- 
perfect layer, which corresponds to that in the k,j.ky mo- 
mentum space, is denoted by . The Green function 
in the system with the imperfect layer and that in the 
bulk system without the imperfect layer are defined as 
G^ and G^, respectively. We approximate the strong 
electron correlation by the on-site repulsion between f- 
electrons in the following Hubbard model: 



H 



(3) 



where c[^{cja) is the creation (annihilation) operator for 
the electron with the spin a at the site i (j). The transfer 
integral denotes the nearest-neighbor hopping ampli- 
tude between the i-th and j-th sites. We set the value of 
the on-site Coulomb repulsion as U — 5.5. We obtain the 
effective interaction for the d-wave pairing on the basis 
of the second-order perturbation theory with respect to 



U. 

In this study, we are interested in the imperfect-layer 
superconductivity, which is bound to the imperfect layer 
owing to the site potential W. The localized states, which 
are called imperfect-layer states, are energetically separa- 
ble from the bulk states and appear under the condition 
that \W\ is larger than ^tz- This condition is obtained 
by the analysis of the Green function. -'^^ When W > ^t^, 
the wave functions of superconductivity, as well as the 
Green function, are evanescent along the z-direction ow- 
ing to the effect of W. Thus, the decay length of the 
Green function is sufficiently short along the z-direction 
in the energy range of predominant excitations of the 
superconductivity. This situation indicates that the su- 
perconducting equations can be separable between the 
imperfect layer and the bulk layers, and then, Tcimp is 
also independent of Tc^buik- Therefore, the effective pair- 
ing interaction between layers can be neglected as an ap- 
proximation, and the superconducting equations for each 
layer are obtained as 



k' ,u 



Gm.m(fc|h'^OG'm,m(-fc||, —Uj'i)Ar 



(4) 



''\\T^U^m,m\~'^\\T-^lJ^m,m{k\\,Uj'i). 

Here, Am indicates an eigenvalue in the TO-th layer. The 
anomalous self-energy Am,m has the (ix2_y2-wave pair- 
ing symmetry because the dx^-y^'-wave pairing is domi- 
nant in the bulk superconductivity of Ce-115 compounds. 
This result has been explained by experimental and the- 
oretical studies. The effective interaction for the 
d-wave pairing in the TO-th layer is given as 

Vrn,mik\\ " k'p^l " ^l) = U + U'^Xm.rnik\\ 

The susceptibilities Xm,n are expressed as 
T 



k',UJl-L0'i). (5) 



Gm,„(fc|| +g||,CJi + f7;)G^_„(fc||,W;)- (6) 



The diagonal component Xm,m corresponds to the sus- 
(fl^ptibilities in the m-th layer. The boson Matsubara fre- 
quency is denoted by fli = 2'kTI and the number of sites 
in the plane is defined as A^|| . The electron state in the m- 
th layer becomes the superconducting state at the tran- 
sition temperature when the eigenvalue Xm of the TO-th 
layer reaches 1. The behavior of Am corresponds to the 
superconducting transition temperature, and thus, the 
large Am corresponds to the high transition temperature. 
We calculate numerically the dependence of Am on to at 
T = 0.0005. In our numerical calculations, we divide the 
first Brillouin zone of k^ky into 128^ meshes and use 2048 
Matsubara frequencies. 

In Fig. 2, we show the dependence of Am on to at T = 
0.0005 as the numerical result. The eigenvalue Aq of the 
imperfect layer is markedly increased in the case that the 
value of the site potential is close to = —1.5tz, whereas 
Am is not increased in the imperfect layer with W = 
and the other values except for W = —1.5tz- The values 
of Am of the other layers do not increase with W, and this 
value is almost equal to that of the bulk eigenvalue. This 
result indicates that Tc^imp becomes higher than Tc_buik 
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Fig. 2. Dependence of A„ 
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Fig. 3. Effect of W on po(<^) in the imperfect-layer state at T : 
0.0005. 
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To clarify the effect of W on Tc^mp, we show the den- 
sity of states po{i^) in Fig. 3. The sharp peak of po{i^) 
locates near the Fermi level a,t W = —1.5tz, and then, it 
induces Tc^imp to be higher than Tc^buik- The sharp peak 
is not present near the Fermi level in the case of = 
and is not close to the Fermi level at the other values ex- 
cept for W = — 1.5iz. As —W becomes larger, the peak 
shifts to the lower energy side and grows sharply. 

We show the Fermi surface of the imperfect-layer state 
in the kxky momentum space in Fig. 4(a). It has a 47.1% 
area occupied by the electrons in the momentum space, 
and thus, the occupation of the Fermi surface corre- 
sponds to the near half-filling involving a strong nesting 
feature. The nesting feature induces the strong antiferro- 
magnetic spin fluctuations at (tt, tt), which appear in the 
susceptibility xo,oi^i — 0) shown in Fig. 5. The momen- 
tum dependence of the antiferromagnetic spin fluctua- 
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Fig. 4. 

momentum space at T = 0.0005. It has the strong nesting fea- 
ture, (b) Fermi surfaces in the fc^^fcy momentum space at = 0, 
7r/2, and tt at T = 0.0005 in the case of the bulk system without 
the imperfect layer. 



tions is in favor of the (ia.2_y2-wave pairing. The above- 
mentioned features of both the strong antiferromagnetic 
spin fluctuations and the large po{^) near the Fermi level 
markedly increase the Am of the imperfect layer, and 
then, the increase in Am leads to that in Tc^mp- The 
features of the imperfect-layer state are similar to the 
d-wave superconductivity in the high-Tc copper oxides^^ 
near half-filling. 

On the other hand, the features of the bulk states are 
different from that of the imperfect-layer state. We show 



the Fermi surfaces of the bulk state in the k^ky momen- 
tum space at kz = 0, tt/2, and tt in Fig. 4(b). The electron 
filling in the bulk state is equal to 0.384. Although the 
Fermi surface at kz = tt has the strong nesting feature 
similar to that of the imperfect-layer state, the nesting 
features are not so strong on the Fermi surfaces at the 
other values of kz including kz = and tt/2. While the 
strong nesting feature at fc^ = tt produces the strong 
antiferromagnetic spin fluctuations leading to the high 
transition temperature, the Fermi surfaces at kz = and 
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Fig. 5. Susceptibilities XO,o('^i = 0) iii the imperfect-layer state 
and the bulk state in the kxky momentum space at T = 0.0005. 
Here, we show the susceptibility Xo.o('»^l = 0) at fcz = in the 
bulk system without the imperfect layer. The antiferromagnetic 
spin fluctuations in the imperfect-layer state are enhanced at 
(tt, 7r) owing to the effect of W . 



7^/2 do not. Therefore, Tc^buik is lower than Tc.imp because 
it is determined by the effective interaction involving the 
entire set of Fermi surfaces at all values of fc^ . 

In our conclusion, the effect of the site potential in the 
imperfect layer induces to be Tc^mp higher than Tc.buik- 
The effect of the site potential produces the electron state 
bound to the imperfect layer, which is independent of 
the bulk states. The Fermi surface of the imperfect-layer 
state has the 2-dimensional nesting feature, which leads 
to increases in both the antiferromagnetic spin fluctua- 
tion and the density of states near the Fermi level, which 
are in favor of the (ia;2_j,2-wave superconductivity. We 
determined the value of the site potential leading to this 
situation and clarified that Tc^mp becomes higher than 
7c, bulk owing to the effect of the site potential in the im- 
perfect layer. From the results, we propose the existence 
of an imperfect layer in Celrlns, which gives T^^p higher 
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